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Abstract: 

We calculate the rate for the decay B u — > lwy using the light cone QCD sume rules. 
We find Br{B u — > lisj) ~ 2 • 10~ 6 . The results are used to test the applicability of 
the constituent quark model approximation to the same process. The latter estimate 
is proportional to l/m^, where m u ~ A u is the "constituent quark mass", indicating 
that the process is of long distance type. We find that the two approaches yield 
similar results for the total rate with the choice m u ~ 480 MeV. This indicates 
that the constituent quark model may be used for estimates of the radiative "an- 
nihilation" contribution to this and other radiative decays. We point out that this 
decay may be useful for the measurement of \V u b\. 
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1 Introduction 



Recently there has been increasing interest in long distance contributions to inclu- 
sive fl], |], |j and exclusive [§], |, §, 0, H radiative D- and B-decays. One way of 
estimating the long distance part of annihilation-type contributions is to use a sim- 
ple constituent quark model in which the photon (or gluon) is emitted from the light 
quark in the initial state |L], |2|, |], |6|, [7|, |8[ . Thus for radiative decays such as B — > /z/7 
or B — > p7 one finds that the annihilation amplitude is proportional to l/m q , where 
q is the light quark in the decaying meson. This indicates that it is indeed of the 
long distance type in an average sense, i.e. it is equivalent to a sum of all possible 
hadronic intermediate states p[. One expects that more detailed observables, such 
as the photon spectrum in B — ► /z/7, will not be reproduced in full detail by the 
quark model. 



In a recent paper [[7] the reaction B — > Iw-f was discussed using the annihilation 
static quark diagram with the photon emitted from the u-quark giving the largest 
contribution. The result of this calculation exhibits the particular dependence 
(see Eq. (25) below). One may expect that in the case at hand this constituent quark 
mass m u must be close to the " inertia parameter" A n of the Heavy Quark Effective 
Theory, which is defined as A = lim mQ ^ 00 (M(Q g ) — thq). Estimates of A n typically 
yield a somewhat larger value for this quantity than the number m u ~ 350 MeV, 
which holds for light hadrons (see e.g. ||). 



In this letter we suggest a different and largely model-independent method to 
calculate the differential and total width of the decay B — > Iwy. Our approach is 
based on the use of the concept of duality in a form suggested by the QCD sum 
rules method [jTIJ. In this approach an amplitude of interest is extracted from the 
imaginary part of a suitably defined Euclidean correlation function. More specifi- 
cally, we use here a version of the QCD sum rule technique known as the light cone 
QCD sum rules. The method has been developed for light quark systems in Ref. 
||1 1||, ||1 2j|. It has proved to be a convenient tool for a study of exclusive processes 
with the emission of a light particle. An Euclidean correlation function correspond- 
ing to such a process can be calculated using the light final particle as a source of a 
properly defined variable external field in which the correlation function develops. 
It turns out that in this case different operators contribute according to their twist 
rather than their dimension. Matrix elements of non-local operators in the variable 
external field are identified with the set of wave functions of increasing twist, and 
replace vacuum expectation values of local operators which appear in the traditional 
sum rules method. The behavior of wave functions is severely restricted by the (ap- 
proximate) conformal invariance of QCD. More details on the method and a list of 
references can be found in [p~3], WM, see also Sect. 2 and 3 below. 



In our problem we deal with a set of the photon wave functions of increasing twist. 
Photon matrix elements are defined via a sum of matrix elements containing vector 



1 



meson states (see sect. 3 for more details). Unlike perhaps more complicated hadron 
(such as 7T- or p-) wave functions, which are still debatable in the literature, the 
photon wave functions turn out to be rather simple and are given by their asymptotic 
expressions. Thus theoretical uncertainties are minimal in this case. We note in 
passing that our approach in not equivalent to the direct use of the vector meson 
dominance applied in || to compare with the quark model annihilation contribution 
in the decay B — > pj. The paper is organized as follows. In sect. 2 we obtain the light 
cone sum rule for the process of interest to the twist 2 accuracy (with partial account 
for twist 4 effects). In sect. 3 we collect the necessary information on the photon wave 
functions of twist 2. The final sum rule and a numerical analysis are presented in 
sect. 4. We find transition form factors and calculate the differential and total decay 
width. The results are then confronted with those of Ref. M. We find that the rates 
obtained by these two methods are the same for the values of the constituent quark 
mass m u ~ 480 MeV. This range is in reasonable agreement with independent 
estimates for A. These results indicate that the quark model approximation may 
be used for estimates of the radiative "annihilation" contribution to this and other 
heavy meson decays, in which the photon is dominantly emitted from the light 
quark, such as B° — > D°*7 , i? ± — > p 7 , B s — > 77. Sect. 5 contains a summary 
and several concluding remarks. We finally comment on the possible usefulness of 
this decay for the measurement of \V u b\ and discuss some general problems related 
to the annihilation mechanism in exclusive amplitudes. 



2 Light cone sum rule 



The effective Hamiltonian for the decay of interest is 

Heff = %V ub ^(1 - 75)/ (7|«7m(1 - 7b)6|5) (1) 

We have neglected here the helicity-suppressed photon emission from the final lepton. 
We use the following parametrization for the hadron transition form factors 

-^=( 7 (g)|ur M 6|B(p)> = e^ pa el Pp q a ^- (2) 



+ i{el{pq)-{e*p)q,) f{(r] 



m 



B 



(here r M = 7^(1 — 75) , Q = p — q and e* stands for the polarization vector of the 
photon.) Our aim is to calculate the transition form factors g , / including their 
momentum dependence. To this end we use the light cone QCD sum rules method. 
Technically, our method is very close to the calculation of the semileptonic B — ► nev 



T5|j and radiative B — > p + 7, B — > K* + 7 [16| form factors. These calculations are 



based on using the light cone wave functions of a light pseudoscalar ( tt) or vector 
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( p, K*) meson, respectively. In our problem, the relevant objects are the light cone 
photon wave functions, which will be discussed below. 



We start with the correlation function 



T,(p,q) = ij <?xe^(l(q)\T{ii(xft f (l-nMx)HO)iw(0)}\<>) 

= ^^W (7(,)|fio>(l - 75)6|B(p+,)>+ --- (3) 

= yfA-na \e fivpa e f v 'p p q (J T x + i(e*(pq) - (e*p)q fl )T 2 } 

where (B\bi~f 5 u\0) = vn^j^l^b . The second line in the above formula selects the 
pole contribution due to the B-meson pole in the corresponding dispersion integrals 
for the scalar form factors g and /. On the other hand, the correlation function (3) 
can be calculated in QCD at large Euclidean momenta (p + q) 2 . 

In this case the virtuality of the heavy quark, which is of order m 2 — (p + q) 2 , 
is the large quantity. Thus, we can expand the heavy quark propagator in powers 
of slowly varying fields residing in the photon, which act as external fields on the 
propagating heavy quark. The expansion in powers of an external field is also the 
expansion of the propagator in powers of a deviation from the light cone x 2 ~ 0. 
The leading contribution is obtained by using the free heavy quark propagator in 
the correlation function (3). We obtain 

(27r)4K _ k2 f {p - k)x (i(q)Wx)^(i - 76 )(vr„ + m b )z 75 n(o)|o) (4) 

In this formula a path-ordered gauge factor between the quark fields is implied, as 
required by gauge invariance. In the particular case of the Fock-Schwinger gauge 
x^A^ix) = this factor is equal to unity. 

We see that the answer is expressed via the one-photon matrix element of the 
gauge invariant non-local operator with a light-like separation x 2 ~ 0. This matrix 
element defines a light cone photon wave function (WF) in a way analogous to the 



definition of the light cone meson WF's fll7| . We define two-particle photon WF's 
as follows : 



(j(q)\u(x)a fw u(0)\0) = iV4jme u (uu) du e luqx [(e^q v - e I/ g M )(x0(w) 

J 

+ x 2 g { - 1 \u)) + {{qx)(e li x„ - e u x^) + (ex)(x fl q v - x u q^) 
- x 2 (e fl q v - e v q^)}g {2 \u)} (5) 



(7(g)|M(x)7 At 7 5 M(0)|0) = -\P&WLe iiUpa e u qpX (J f / due luqx g±(u) (6) 

4 Jo 

Here <j){u) , g±(u) stand for the leading twist 2 photon WF, while and g^ are the 
two-particle WF's of twist 4 (see below for more detail). The above WF's describe 
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a distribution in the fraction of the total momentum carried by the quark (up z ) and 
the anti-quark ((l — u)p z ) in the infinite-momentum frame p z — > oo (or, equivalently, 
a distribution in the light cone momentum P+). The dimensional constants x an d 
/ are chosen in such a way that a function f{u) = {4>{u) , g±(u)) is normalized to 
unity : 

f 1 duf(u) = l (7) 
Jo 

It is assumed that perturbative logs of x 2 are resummed by renormalization group 
methods and result in a dependence of the WF's on the scale factor u 2 ~ x~ 2 ~ 
m b ~ (P + ( l) 2 i w hich in other words fixes the normalization point of local operators 
arising in an expansion of the non-local matrix elements (5-7). The corresponding 
anomalous dimensions can be found in the literature. 



The WF's g( ' , represent twist 4 contributions to the two-particle photon 
WF (5). Using equations of motion one can relate them with three paricle WF's of 
twist 4 which include an additional gluon [|ll|, 12| . These three-particle WF's would 
appear in the correlation function (3) if a soft gluon emission from the heavy quark is 
taken into account. Such a contribution is typically rather small and will be omitted 
in what follows. On the other hand, we retain the two-particle twist 4 contribution 
due to g^ l > , g( 2 \ That is, the most important twist 4 effects are included in our 
calculation. Finally, adding a perturbative photon emission from the u- and b- 
quarks we arrive at the following expressions for the invariant functions T\ , T 2 : 



To 



+ 



+ 



1 



o mf — (p + uq) 2 



g±{u) 



— f 

2 ml — (p + uq) 2 
1 1 
o ml — (p + uq) 2 



e u (uu) \ X<P{u 
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We now proceed to a description of the photon WF's (5-6). 



3 Light cone photon wave functions 



The leading twist photon WF (f>(u, u 2 ) has been introduced and investigated in detail 
in Ref.(Tl|] in connection with a study of the radiative decay S — > p+'j. Formally, the 
photon matrix element of an arbitrary (non-local) operator is defined as a weighted 



sum of matrix elements of the same operator between the vacuum and (transversely 
polarized) vector meson states : 

(7(?)|...|0> = ^v / 4^e( A ) J dz e^'(0\Tj^(z) . . .\0) 

= vW£-(V n (g)|...|0>, (10) 

n 9n 

where j p m = ^uj^u — ^d r y^d — |s7^s and the couplings g n are defined via 

m m \ V n) = i X) — . ^5.5 (11) 
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The particular value for g p quoted above has been found long ago by the QCD sum 
rules method HT0|. 



Our definition (10) is equivalent to the one given in fLl |, [ 12] for the case of 
the leading twist photon WF <p±_(u). There it has been defined as the vacuum 
expectation value of the non-local operator qa^q (where q = u, d, s ) in an ex- 
ternal electromagnetic field F a p(x). The QCD action is modified in this case to 
/ dz Lqcd(z) + V^na J dz A fl (z)j^ n (z) and, expanding to the first order in the plane 

wave field strength F a p(x) = i(epq a — e^q^)e iqx , we obtain 



dz e^(0\T{f p m (z)q(x)a a(3 q(0)}\0) = e q (qq) X (qpS aft - q a S P(t ) J du e^0(«) (12) 
Thus, according to Eq.(10), one gets 

(7(g)|M(x)a Q/3 n(0)|0) = i^^e u {uu) X {e^qp - efq a ) f du e^^u) , (13) 



i.e. Eq.(5). In the constant external field limit q — > 0, Eq.(12) goes back to the 
definition of a quantity known as the magnetic susceptibility of the quark condensate 

(qcr a 0q)F = e q (qq)xF a p (14) 
The magnetic susceptibility x nas been found pO] , pl| by the QCD sum rules 



method. At an accuracy of the order of 30 % it is dominated by the p-meson 
contribution : at p? ~ 1 GeV 2 

_ J —4.4 Gev~ 2 p, p', p" are included 
1 —3.3 Gev~ 2 p is included 

It has been suggested on these grounds that the vector dominance works sufficiently 
well for electromagnetic properties of hadrons in the long wave limit q — > 0. Thus, to 
the quoted accuracy, by virtue of the local version of Eq.(10) we obtain the following 
value for the normalization constant / in Eq. (6) : 

f~^f p m p , f p ~ 200 MeV (15) 
9p 



The validity of Eq.(15) can be checked by the standard QCD sum rules technique. 



Here some comments are in order. It is known that two-particle operators con- 
taining 7^ and 7^75 instead of o^ v have non-zero projections on a vector meson 
state. The corresponding wave functions ( call them g\ and g±, respectively) con- 
tain pieces of different twist. As has been shown in Ref. [lfj, a twist 2 contribution 



to a vector meson WF g\ can be expressed via an integral of a longitudinally po- 
larized vector meson WF of the leading twist 2. A similar relation holds for the 
derivative dg±(u)/du (see Eq.(39) of ). Thus we conclude that in the case of 
photon the WF g v ± vanishes to the leading twist accuracy, while the WF g\ is a con- 
stant which is set to be unity with our normalization. At this point we differ from 



similar calculations in Refs. |23|, ^4J where a contribution due to g± was omitted. 
Numerically, the effect of the inclusion of the corresponding term in the sum rule 
(see Eq. (19) below) is about 25 % for the decay width. 

Let us now discuss the much less trivial problem of finding the functional in- 
dependence of the WF's (5-6). As a non-local operator is equivalent to an infinite 
series of local operators , this task might appear formidable. Indeed, in order to re- 
store the functional w- dependence, one has to know all moments Jq 1 du (2u — l) n f(u) 1 
i.e. , according to (5-6), all matrix elements of the type ( r y\q( y \/ a x a ) n a^ u q\0) . The 
way out lies in the use of an approximate conformal symmetry of QCD which holds 
at the one-loop level. The conformal invariance permits an expansion of a non-local 
operator of the type (5-6) in a series over multiplicatively renormalizable operators 
with ordered anomalous dimensions, which induces the expansion of the leading 
twist WF in the series of Gegenbauer polynomials ||i~l|,[12"f,|18|. The asymptotic 



WF's are defined as contributions of operators with the lowest conformal spin and 
unambiguously fixed by the group structure. Pre-asymptotic corrections correspond 
to the operators with the next-to-leading conformal spin, whose numerical values 



are usually calculated by the QCD sum rules method. It was found |L1[ that for 
the leading twist photon WF contributions of higher conformal spins are small. In 
terms of the dispersion relation this means that the corrections to the asymptotic 
WF's of the p, p', p" states have opposite signs and nearly cancel in the sum. Thus, 
the simple asymptotic formula holds (in what follows, we use u = 1 — u) 

<f>±(u) = 6uu (16) 

To the leading twist accuracy we use for the WF (6) g±(u) = 1 . For the twist-4 
WF's g^(u) , g^ 2 \u) we use the following expressions (see |2c 



«(«) = --u(3-u) 



9 y \ u 

g (2 \u) = (17) 
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4 Form factors and the decay width 



To match the answers (8,9) with the B-meson contribution to the correlation func- 
tion (3), we note that (8,9) can be re-written as the dispersion integrals with the 
expression {m 2 — up 2 )/u being the mass of the intermediate state. The duality pre- 
scription requires that this invariant mass has to be restricted from above by the 
duality threshold s — 35 GeV 2 (this value is obtained from corresponding two-point 
sum rules). As it is easy to see, this transforms into an effective cut-off in the lower 
limit of the u-integral [I5],[I3||. Finally, we make the standard Borel transformation 
suppressing both higher states resonances and higher Fock states in the full photon 
wave function. Under the Borel transformation — (p + q) 2 — > M 2 



1 



m% — (p + q) 2 
1 1 



exp 



m l ~ (P + uc t) 2 
Our final sum rules take the form 
777.6 f 1 du ,m° 
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e u (uu) ( X <t>(u) -4(p«(«) -g^{u))^ 



ml + uM 2 ' 
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(e u - e b )u- 
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mt 



+ e b In • 
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f B Jo u M 2 
e u (uu) x0(m) -Ag [l) {u) 



uM 2 
ml + uM 2 
u 2 M A 



s -p* 




e u - e b ){u - u + 



p- 



p 2 u 2 



+ 



m b f 
2uM 2 



(19) 



rnr 



mt 



up 2 ' 



\ p 2 \ u(m 2 b -p 2 ) _ p 2 m 2 b - up 2 
(e u + e b )u— — 2 ~1T + e b\ u ~ u H 2) ln ' 



m 6 / m 6 — M P 



mt 



umt 



(20) 



We now turn to numerical estimates. In evaluating (19,20), we have used the fol- 



lowing set of parameters 
135 MeV [I 



m b 



4.7 GeV, m B = 5.28 GeV, s ~ 35 GeV 2 , f. 



B 



The particular value of fs corresponds to the result of the sum 
rule calculation with 0(a s ) corrections omitted, cf. jn|, |1^]. Within our accuracy 
we also neglect the logarithmic evolution of the WF's to a higher normalization 
point in the sum rules (19,20) which is of order of the characteristic virtuality of the 
heavy quark in the B-meson. The Borel mass M 2 has been varied in the interval 
from 8 to 20 GeV 2 . We have found that within the variation of M 2 in this region, 
the result changes by less than 10 %. As the sum rules of the type (19,20) are 
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expected to work in the region m\ — p 2 ~ a few GeV 2 , which is somewhat smaller 



than the maximal available p 2 = m 2 , we have to use some extrapolation formulas 



to extent the results of the sum rules to the whole region of p 2 . We have found that 
the best agreement is reached with the dipole formulas 

g(p 2 )^ ^— , f(p 2 )^ ^— (21) 

\ rrvj ' \ m| / 

The fact that the dipole approximation agrees better with results of the sum rules at 
not too large p 2 than the pole-type formulas has been also noted in [16| for analogous 
sum rules including vector mesons. We find 

h ~ 1 GeV , mi ~ 5.6 GeV 
h 2 ~ 0.8 GeV , m 2 ~ 6.5 GeV (22) 



The calculation of the decay width yields 

r 



aG 2 F \V ub \ 2 m% 

96vr 2 1 ' 



where (x — 1 — x) 



mil = f 1 dxx 3 x{g 2 (x) + f{x)} (24) 
Jo 

Here x = 1 — 2E/ms and i£ stands for the photon energy. The corresponding 
differential width can be read off Eqs. (23,24). We expect the accuracy for / to be 
of the order of 20-30 %. The numerical answer for the quantity / is I ~ 0.021 
( I ~ 0.015 if the contribution due to g± is omitted). 

We would like now to compare our result with that obtained in Ref. within 
the non-relativistic quark model. These authors have obtained the following answer 
for the decay width : 

_ aG 2 F \V ub \ 2 m% Q 2 U f 2 
487r^ 6 m„ 

where now fs stands for the physical value fs — 175 MeV and a small contribution 
stemming from photon emission off the b-quark has been neglected. Note that in 
Eq. (23) the constant fs appears implicitly in the denominator (see Eqs. (19,20)), 
while it is in the numerator of Eq.(25). This difference stems from the use of the 
dispersion approach. We find that Eq.(25) yields a result numerically close to (23) 
with the choice of the free parameter 

m u ~ 480 MeV , (26) 

The photon spectra corresponding to the light cone QCD calculation and the quark 
model approach are compared in Fig.l using the same normalization. One sees 
that while the quark model gives a fully symmetric spectrum, in our approach the 
spectrum is slightly asymmetric, as a result of a balance between a typical highly 
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asymmetric resonance-type behaviour given by the non-perturbative contributions, 
and a perturbative photon emission. Both curves emphasize the dominant role of 
hard photons in the decay of interest and vanish at the end-points. 



Finally, we present the estimated rate for the decay of interest 

Br(B -> J1/7) ~ 2 ■ 10~ 6 , (27) 
which corresponds to the choice m u = 480 MeV. We have used here the values 



\V ub \/\V cb \ = 0.08 , \Vctl = 0.04 |2 



5 Summary 



We have calculated within the light cone QCD sum rules technique the transition 
form factors and the differential and total width of the radiative decay B — > lwy. 
In accordance with the expectations of Ref. j7|, the photon emission overcomes the 
helicity suppression of the pure leptonic B — > Iv amplitude and yields an exper- 
imentally admissible decay width. Comparing our result with that suggested by 
the quark model, we have fixed the value of the effective constituent quark mass 
m u ~ 480 MeV, which is consistent with estimates for the "inertia parameter" A u 
in the heavy quark effective theory. The photon spectra calculated within our ap- 
proach and the quark model calculation have similar shapes. They both emphasize 
hard photon emission, vanish at the boundaries of the spectrum and do not exhibit 
the typical bremsstrahlung 1/E behavior. The completely symmetric form given 
by the quark model is an artifact of the absence of hadronic poles. The light 
cone QCD approach gives a slightly asymmetric spectrum due to a balance between 
resonance-type non-perturbative contributions and perturbative photon emission. 

In the perturbative QCD approach to exclusive processes [JT7J annihilation am- 
plitudes are typically expressed via a convolution of wave functions of participating 
hadrons. These integrals turn out to be logarithmically divergent in the end-point 
region |T7|], which indicates the breakdown of perturbative factorization and sug- 
gests that the annihilation mechanism is of long distance type. This phenomenon 
is a manifestation of a long distance dominance, akin to an appearence of the l/m q 
factor in the quark model approach. In this sense the situation is similar to the 



problem of the Feynman-like end-point contribution to the pion form factor |]T3 
We believe that the dispersion light cone QCD approach, suggested in this context 
in Ref. [O and used in the present letter for the particular annihilation-like process, 



may be of use also for reactions involving gluon annihilation. 

We conclude by remarking that this process may be useful to extact the KM 
matrix element \V u b\ , since the value obtained for the branching ratio Br(B — ► 
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fii"y) ~ 2 ■ 10 6 is reasonably large and, in addition, is theoretically under control 
in view of the simplicity of the hadronic system. 

Note added 

After completion of the previous version of this work we became aware of recent 
papers |23], [24[] where a similar technique has been applied to the radiative decays 
B — > p + 7 and B — > uj + 7 and in |24]] to the B — > /z/7 decay. 



The research of G.E. was supported in part by GIF and by the Fund for the 
Promotion of Research at the Technion. The work of R.R.M. was supported in 
part by the Natural Sciences and Engineering Research Council of Canada. We are 
grateful to V.M. Braun for pointing out errors in an earlier version of this work 
where we used an incorrect set of photon wave functions of leading twist and did 
not take into account perturbative photon emission. 
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Figure 1: The normalized differential width A = l/T(dT/dx) of the decay B — > Ivy 
as functions of x — 1 — 2E/m,B where E is the photon energy. The solid and dashed 
curves correspond to the results of the light cone sum rule calculation (Eqs. (25,26)) 
with the dipole ansatz for the form factor f(p 2 ) and the quark model calculation of 
Ref. J7|, respectively 



13 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9506264v2 



